29 research outputs found
AMBER: A Semi-Numerical Abundance Matching Box for the Epoch of Reionization
The Abundance Matching Box for the Epoch of Reionization (AMBER) is a
semi-numerical code for modeling the cosmic dawn. The new algorithm is not
based on the excursion set formalism, but takes the novel approach of
calculating the reionization-redshift field
assuming that hydrogen gas encountering higher radiation intensity are
photoionized earlier. Redshift values are assigned while matching the abundance
of ionized mass according to a given mass-weighted ionization fraction
. The code has the unique advantage of allowing users to
directly specify the reionization history through the redshift midpoint
, duration , and asymmetry
input parameters. The reionization process is further controlled through the
minimum halo mass for galaxy formation and the radiation mean
free path for radiative transfer. We implement improved
methods for constructing density, velocity, halo, and radiation fields, which
are essential components for modeling reionization observables. We compare
AMBER with two other semi-numerical methods and find that our code more
accurately reproduces the results from radiation-hydrodynamic simulations. The
parallelized code is over four orders of magnitude faster than radiative
transfer simulations and will efficiently enable large-volume models, full-sky
mock observations, and parameter-space studies. AMBER will be made publicly
available to facilitate and transform studies of the EoR.Comment: 29 pages, 21 figures, 1 table. Submitted to ApJ. AMBER will be made
publicly available when the paper is publishe
Triple and Quadruple Black Holes in the ASTRID Simulation at
We use the ASTRID cosmological hydrodynamic simulation to investigate the
properties and evolution of triple and quadruple Massive Black Hole (MBH)
systems at . Only a handful of MBH tuple systems have been detected to
date. In ASTRID, we find of the are in tuples
with . The tuple systems span a range of
separations with the majority of the observable AGN systems at kpc. They include some of the most massive BHs (up to ) but with at least one of the components of . Tuples' host galaxies are typically massive with . We find that massive halos with host MBH tuples. Following the subsequent interactions between
MBHs in tuples, we found that in of the triplets all three MBHs
merge within a Gyr, and go through one merger. As a by-product of the
complex multi-galaxy interaction of these systems, we also find that up to
of tuples lead to runaway MBHs. In ASTRID, virtually all of the
ultramassive black holes () have undergone a triple quasar
phase while for BHs with this fraction drops
to .Comment: 10 pages, 9 figures; comments welcom
Fly-by galaxy encounters with multiple black holes produce star-forming linear wakes
We look for simulated star-forming linear wakes such as the one recently
discovered by van Dokkum et al. (2023) in the cosmological hydrodynamical
simulation ASTRID. Amongst the runaway black holes in ASTRID, none are able to
produce clear star-forming wakes. Meanwhile, fly-by encounters, typically
involving a compact galaxy (with a central black hole) and a star-forming
galaxy (with a duo of black holes) reproduce remarkably well many of the key
properties (its length and linearity; recent star formation, etc.) of the
observed star-forming linear feature. We predict the feature to persist for
approximately 100 Myr in such a system and hence constitute a rare event. The
feature contains a partly stripped galaxy (with ) and a dual BH system () in
its brightest knot. X-ray emission from AGN in the knot should be detectable in
such systems. After from the first fly-by, the
galaxies merge leaving behind a triple black hole system in a (still) actively
star-forming early-type remnant of mass .
Follow-up JWST observations may be key for revealing the nature of these linear
features by potentially detecting the older stellar populations constituting
the bright knot. Confirmation of such detections may therefore help
discriminate a fly-by encounter from a massive BH wake to reveal the origin of
such features.Comment: 8 pages, 5 figures, comments welcom
Orbital and Radiative Properties of Wandering Intermediate-Mass Black Holes in the ASTRID Simulation
Intermediate-Mass Black Holes (IMBHs) of are commonly
found at the center of dwarf galaxies. Simulations and observations
convincingly show that a sizable population of IMBHs could wander off-center in
galaxies. We use the cosmological simulation ASTRID to study the orbital and
radiative properties of wandering IMBHs in massive galaxies at . We
find that this population of black holes has large orbital inclinations
() with respect to the principal plane of the host. The
eccentricity of their orbits is also significant () and decreases
with time. Wandering IMBHs undergo spikes of accretion activity around the
pericenter of their orbits, with rates times the Eddington
rate and a median accretion duty cycle of . Their typical spectral
energy distribution peaks in the infrared at rest-frame.
Assuming a standard value of for the matter-to-energy radiative
efficiency, IMBHs reach keV X-ray luminosities for of the time. This luminosity corresponds
to fluxes within Mpc.
They could be challenging to detect because of competing emissions from X-ray
binaries and the interstellar medium. X-ray luminosities , in the hyper-luminous X-ray sources (HLXs) regime, are
reached by of the IMBHs. These findings suggest that HLXs are a
small subset of the wandering IMBH population, which is characterized by
luminosities times fainter. Dedicated surveys are needed to assess
the demographics of this missing population of black holes.Comment: Accepted for publication in MNRAS. This is the final version of the
manuscript. 9 pages, 7 figure
A vast population of wandering and merging IMBHs at cosmic noon
Massive black holes in the centers of galaxies today must have grown by
several orders of magnitude from seed black holes formed at early times.
Detecting a population of intermediate mass black holes (IMBHs) can provide
constraints on these elusive BH seeds. Here we use the large volume,
cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and
dynamical friction to investigate the population of IMBH seeds. Dynamical
friction is largely inefficient at sinking and merging seed IMBHs at high-z.
This leads to an extensive population (several hundred per galaxy) of wandering
IMBHs in large halos at z~2. A small fraction of these IMBHs are detectable as
HLXs, Hyper Luminous X-ray sources. Importantly, at z ~ 2, IMBHs mergers
produce the peak of GW events. We find close to a million GW events in Astrid
between z=2-3 involving seed IMBH mergers. These GW events (almost all
detectable by LISA) at cosmic noon should provide strong constraints on IMBH
seed models and their formation mechanisms. At the center of massive galaxies,
where the number of IMBHs can be as high as 10-100, SMBH-IMBH pairs can form.
These Intermediate mass ratio inspirals (IMRIs) and extreme mass ratio
inspirals (EMRIs), will require the next generation of milli-muHz space-based
GW interferometers to be detected. Large populations of IMBHs around massive
black holes will probe their environments and MBH causal structure
A Close Quasar Pair in a Disk-Disk Galaxy Merger at z = 2.17
Most local massive galaxies, if not all, are believed to harbor a
supermassive black hole (SMBH) at the center. Galaxy mergers have long been
thought to drive strong gas inflows and accretion onto one or both central
SMBH, triggering single or dual quasars as a natural stage of the hierarchical
galaxy and SMBH evolution. While many dual active galactic nuclei -- the
low-luminosity counterparts of quasars -- have been observed at low redshift,
no unambiguous dual quasar is known at cosmic noon (z>~2) when both quasar
activity and global star formation density peaked. While a handful of dual
quasar candidates were known at z>1, competing explanations remained. Here we
report multi-wavelength observations of SDSS J0749+2255 as the first kpc-scale
dual quasar confirmed to be hosted by a galaxy merger at cosmic noon. Hubble
Space Telescope NIR imaging reveals extended host galaxies underlying the
compact double nuclei (separated by 0.46" or 3.8 kpc) and tidal features as
evidence for galactic interactions. We also present new multi-wavelength
observations, all lending support to the dual quasar hypothesis. Unlike the
low-redshift low-luminosity counterparts, the high-redshift dual quasar is
hosted by two massive compact disk-dominated galaxies, which may be critical
for efficient gas fueling onto the SMBHs in the early-stage merger. The
apparent lack of stellar bulges and that SDSS J0749+2255 already follows the
local SMBH mass-host stellar mass relation are at odds with the canonical
SMBH-host co-evolution picture and suggest that at least some SMBHs may have
formed before their host stellar bulges. While still at kpc-scale separations
where the host-galaxy gravitational potential dominates, the SMBHs may evolve
into a gravitationally bound binary system in ~0.22 Gyr. The merger products at
low redshift are expected to be gravitational wave sources for pulsar-timing
arrays (abridged).Comment: 79 pages, 17 figures, 6 tables; submitte
Massive Black Hole Mergers with Orbital Information: Predictions from the ASTRID Simulation
We examine massive black hole (MBH) mergers and their associated
gravitational wave signals from the large-volume cosmological simulation
Astrid. Astrid includes galaxy formation and black hole models recently updated
with a MBH seed population between and and a sub-grid dynamical friction (DF) model to follow the MBH
dynamics down to . We calculate initial eccentricities of
MBH orbits directly from the simulation at kpc-scales, and find orbital
eccentricities above for most MBH pairs before the numerical merger.
After approximating unresolved evolution on scales below , we find that the in-simulation DF on large scales accounts
for more than half of the total orbital decay time ()
due to DF. The binary hardening time is an order of magnitude longer than the
DF time, especially for the seed-mass binaries ().
As a result, only of seed MBH pairs merge at after
considering both unresolved DF evolution and binary hardening. These
seed-mass mergers are hosted in a biased population of galaxies with the
highest stellar masses of . With the higher initial
eccentricity prediction from Astrid, we estimate an expected merger rate of
per year from the MBH population. This is a factor of
higher than the prediction using the circular orbit assumption. The LISA events
are expected at a similar rate, and comprise seed-seed mergers,
involving only one seed-mass MBH, and mergers of
non-seed MBHs.Comment: 17 pages, 13 Figures; comments are welcom
PRIYA: A New Suite of Lyman-alpha Forest Simulations for Cosmology
We present the PRIYA suite of cosmological simulations, based on the code and
hydrodynamic model of the ASTRID simulation, and designed for cosmological
analyses of the Lyman- forest. Our simulation suite spans a
-dimensional parameter space, including cosmological parameters and
astrophysical/thermal parameters. We have run low fidelity simulations
with particles in a Mpc/h box and high fidelity simulations
with particles in a Mpc/h box. All our simulations include a
full physics model for galaxy formation, including supernova and AGN feedback,
and thus also contain a realistic population of DLAs. We advance on earlier
simulations suites by larger particle loads, by incorporating new physical
models for patchy hydrogen and helium reionization, and by self-consistently
incorporating a model for AGN feedback. We show that patchy helium reionization
imprints an excess in the 1D flux power spectrum on large scales, which may
allow future measurements of helium reionization bubble sizes. Simulation
parameters are chosen based on a Latin hypercube design and a Gaussian process
is used to interpolate to arbitrary parameter combinations. We build a
multi-fidelity emulator for the 1D flux power spectrum and the mean IGM
temperature. We show that our final interpolation error is and that our
simulations produce a flux power spectrum converged at the percent level for
- . Our simulation suite will be used to interpret Lyman-
forest 1D flux power spectra from SDSS and future DESI data releases.Comment: 24 pages, 11 figures, submitted to JCA
Statistics of Galactic-Scale Quasar Pairs at Cosmic Noon
The statistics of galactic-scale quasar pairs can elucidate our understanding
of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty
cycles of quasar activity in mergers, or even the nature of dark matter, but
have been challenging to measure at cosmic noon, the prime epoch of massive
galaxy and SMBH formation. Here we measure a double quasar fraction of integrated over arcsec separations
(projected physical separations of at ) in
luminous () unobscured quasars at
, using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The
measurement was based on a sample of 60 Gaia-resolved double quasars (out of
487 Gaia pairs dominated by quasar+star superpositions) at these separations,
corrected for pair completeness in Gaia, which we quantify as functions of pair
separation, magnitude of the primary, and magnitude contrast. The double quasar
fraction increases towards smaller separations by a factor of over
these scales. The division between physical quasar pairs and lensed quasars in
our sample is currently unknown, requiring dedicated follow-up observations (in
particular, deep, sub-arcsec-resolution IR imaging for the closest pairs).
Intriguingly, at this point the observed pair statistics are in rough agreement
with theoretical predictions both for the lensed quasar population in mock
catalogs and for dual quasars in cosmological hydrodynamic simulations.
Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST
and Roman, combined with targeted follow-up observations, will conclusively
measure the abundances and host galaxy properties of galactic-scale quasar
pairs, offset AGNs, and sub-arcsec lensed quasars across cosmic time.Comment: 19 pages, 9 figures; submitted to Ap
Tracking SMBH mergers from kpc to sub-pc scales with AXIS
Pairs of active galactic nuclei (AGN) are observational flags of
merger-driven SMBH growth, and represent an observable link between galaxy
mergers and gravitational wave (GW) events. Thus, studying these systems across
their various evolutionary phases can help quantify the role mergers play in
the growth of SMBHs as well as future GW signals expected to be detected by
pulsar timing arrays (PTAs). At the earliest stage, the system can be
classified as a "dual AGN" where the SMBHs are gravitationally unbound and have
typical separations <30 kpc, and at the latest stage the system can be
classified as a "binary AGN" where the two massive host galaxies have likely
been interacting for hundreds of megayears to gigayears. However, detecting and
confirming pairs of AGN is non-trivial, and is complicated by the unique
characteristics of merger-environments. To date, there are less than 50 X-ray
confirmed dual AGN and only 1 strong binary AGN candidate. AXIS will
revolutionize the field of dual AGN: the point-spread-function (PSF),
field-of-view (FOV), and effective area (Aeff) are expected to result in the
detection of hundreds to thousands of new dual AGN across the redshift range 0
< z < 4. The AXIS AGN surveys will result in the first X-ray study that
quantifies the frequency of dual AGN as a function of redshift up to z = 3.5.Comment: 17 pages, 5 figure